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Abstract 

The perturbative QCD approach to multiparticle production assuming Local 
Parton Hadron Duality (LPHD) and some recent results are discussed. Fi- 
nite asymptotic scaling limits are obtained for various observables, after an 
appropriate rescaling, in the Double Logarithmic Approximation (DLA). Non- 
asymptotic corrections are also known in some cases. The DLA applies also 
to very soft particle production where energy conservation constraints can 
be neglected. In this region the particle density follows rather well a scaling 
behaviour over the full energy range explored so far in e + e~ annihilation. 
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1 Introduction 



The study of the intrinsic structure of particle jets produced in hard collisions continues 
to be an active field of research. The interest is directed in elaborating and testing 
specific predictions of perturbative QCD on the parton cascade evolution and secondly, 
to investigate the hadronization process which cannot be treated within a perturbative 
scheme. A reduction of the flexiblity of the models involved and a deeper understanding 
of the phenomenological aspects of the confinement process is an important aim of this 
research. 

The most popular models for particle production in hard collision processes are based 
on a primary hard partonic sub-process which is accompagnied by gluon bremsstrahlung. 
The evolution of the partonic jets is derived in perturbation theory and is terminated at 
a scale of around 1 GeV; thereafter, non-perturbative processes take over and the final 
hadronic particles, often through intermediate resonances, are produced, for example, by 
a string mechanism U or through cluster formation ||. 

Another approach is based on the concept of "Local Parton Hadron Duality" (LPHD) 
0. It has been observed at first that the hadronic energy spectra are rather well rep- 
resented by the parton spectra themselves - without an additional hadronization phase 
- provided the cut-off of the parton cascade is lowered to a value around 250 MeV, of 
the order of the hadronic masses. This general idea has been applied to various other 
observables; the theoretical calculations are based in the simplest case on the Double 
Log Approximation (DLA) |4]||], which provides the high energy limit, or the Modified 
Leading Log Approximation (MLLA) |J which includes finite energy corrections which 
are usually essential to obtain quantitative agreement with experiment at present ener- 
gies 0. Recent experimental results from LEP, HERA and TEVATRON gave further 
support to this approach ||. Although a justification of the model is not yet available at 
a fundamental level the related phenomenology is quite attractive because of its intrinsic 
simplicity with very few parameters. Also the analytical computations allow the deriva- 
tion of scaling laws and the systematics of their violation which provides an important 
insight into the structure of the theory. On the other hand, it is clear that this model 
cannot compete with the standard hadronization models in the description of the various 
details of the final state like the production of different particle species or resonances. It 
has so far been applied successfully for suitably averaged quantities. 

In this presentation we summarize, how scaling and scale breaking predictions obtained 
from analytical calculations compare with experiment. 



2 Basic ingredients of analytical calculations 

We consider high energy collisions which involve a hard subprocess. The colour charges 
of the primarely produced partons are the sources of subsequent gluon bremsstrahlung 
which leads to the partonic jets. The subprocess is described by the corresponding matrix 
element. The gluon bremsstrahlung at small angles 5 with energy E off the primary hard 
parton of type A (A = q, g) with momentum P is given by 
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C A 2 dSdE 2 2N c a s (k T ) f3 2 
dn A =— l0 (k T ) T -, l0 (k T ) = = ^j^ y k T >Q (1) 

where kr ~ ES, (3 2 = 4Nc/b, b = (UNc — 2n/)/3 with Nc, rif the numbers of colours 
and flavours, also C g = N c , C q = 4/3. Inside the cascade the soft gluons are coherently 
produced from all harder partons. For azimuthally averaged quantities the consequences 
of the coherence effect can be taken into account by the angular ordering prescription || 
which requires the angles of subsequent gluon emissions to be in decreasing order. 

The multiparticle properties of the jet can be discussed conveniently by using the gen- 
erating functional fl(| Z A (P, O; u(k)). Here P and denote the initial parton momentum 



and opening angle of the jet, and u(k) is a profile function for particle momentum k. The 
functional is constructed from all the exclusive final states. Then the inclusive densities 
can be obtained by functional differentiation with respect to the profile function u(k) 

p (n \k 1 ,...,k n )=6 n Z{u}/6u(k 1 )...6u(kn) \u=i ■ (2) 

Properties of these densities can be obtained from the evolution equation for Z which 
relates the functional at scales P, O to the one at lower scales according to the "decay" 
A — > BC. In MLLA accuracy this evolution equation is given by || 
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z A ( P , ©) = w r 1 dz 
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(3) 



x ^2 [Zb{zP, 0) Z ((l - z)P, 0) - Z A (P, 0)] 



where ^^(z) denotes the DGLAP splitting functions. The initial condition of the evo- 
lution is given by 

Z A (P,Q;{u})\ Pe = Qo = u A (k = P), (4) 

i.e. at threshold there is only the primary parton. 

These equations take into account energy conservation by choosing the proper argu- 
ments of Z. One can obtain for various observables O analytic solutions which corre- 
spond to the summation of the perturbative series in leading double logarithmic order, 
i.e. the summation of the terms a™L 2n with a large logarithm L. Also results with 
resummed next-to-leading order terms a n s L 2n ~ x are available in some cases. At high 
energies O ~ exp J Y ^y(a s (y)) dy where the anomalous dimension 7 has the expansion 
7 ~ yfa^ + a s + . . . The leading term is refered to as the DLA, the next-to-leading one as 
the MLLA result. The evolution equation (|3]) yields the complete results for the first two 
terms in this ^JoT s expansion. These leading terms are not sufficient, however, to satisfy 
the initial condition (^), this is only possible by using the full result from the summed 
perturbative series. From (0) one can obtain the evolution equations for particle densities 
by appropriate differentiation (^|). Therefore this equation is the basic tool for deriving 
the multiparticle properties of a jet analytically. 

For very high energies the small z contributions dominate (fixed scales Qq and A), 
then one can neglect the recoil effects and approximate 1 — z ~ 1 in the argument of Z 
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in (||). Furthermore, in the high energy limit it is sufficient to include the most singular 
terms of the splitting functions ~ 1/z, in particular, one can neglect the production 
of quark pairs in the cascade with nonsingular splitting function. In this case Eq. (|3]) 
simplifies and can be integrated using the initial condition (|]) to 



Z A (P, 6, it) = u(P) exp ( J dn A [u(E)Z g (E, 6) - 1] 



(5) 



with integration measure from ([I]) and boundary F which takes into account the angular 
ordering constraint 5 < 6 and the cutoff E5 > Q . This is the evolution equation in 
DLA accuracy appropriate for the high energy asymptotics. 



For the energy spectra [jj 1] and a large class of angular correlations []12j one can derive 
from (||) by functional differentiation and appropriate partial integration an evolution 
equation of the type 



h n (S, 9, P) = d n (5, P) + J ^ J ^ 7o K(S, V, K) 



(6) 



where h n denotes generically one such distribution or correlation of order n and d n the 
appropriate initial condition. The singularities in the kernel are regularized by the kx cut- 
off. A nonsingular evolution equation is obtained by changing to logarithmic momentum 
and angular variables. 

3 High energy asymptotics 

The high energy behaviour can be obtained from (|5|). As will be shown in several cases, 
the observable quantities, after appropriate rescaling, approach a finite scaling limit. 



3.1 Multiplicity distribution 



A well known example of such behaviour is the "KNO-scaling" [|i3[|i4| of the multiplicity 
distribution 



< n > PJs) 



tj) = n/ < n(s) > 



(7) 



Here f{ip) is the high energy limit of the probability P n to produce n particles at cms 
energy a/s, rescaled by the average multiplicity < n >. This scaling law has been derived 
in the DLA for the partons of QCD but holds more generally for a large class of 
branching processes [|1^|T^|. Specifically for QCD one obtains also an explicit prediction 
for the function f(ip) or the normalized factorial moments =< n{n — \)...{n — k + 
1) >/< n > k of the multiplicity distribution. For example, for e + e~ annihilation one finds 
F^ = y . This prediction is infrared safe, i.e. independent of the cut-off Qq. Furthermore 
QCD predicts the approach to the asymptotic limit, so for the same quantity one obtains 
with inclusion of the next-to-leading correction (MLLA) HH 



f^ = -a 
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(for rif = 5) which turns out to be large (about 30%); if yet higher order corrections 
are included the result fits the experimental data [18[]. The energy dependence of in 
(0) is very weak (~ 1 / Vhi s) thus simulating the scaling behaviour observed at present 
energies. The ultimate KNO scaling function according to (0) or (|) is broader and will 
be approached only at much higher energies than available today and in the near future. 



3.2 Momentum spectra 

Spectra in the rescaled Bjorken or Feynman momentum variables x = p/P do not scale in 
QCD. Rather the distributions of certain rescaled logarithmic variables approach a finite 
asymptotic limit in the DLA. Such scaling properties have been discussed recently in some 
detail for angular correlations |T2[ (see below). For the energy spectra a scaling limit of 



this type has been suggested already some time ago |ITJ and one obtains in logarithmic 
variables after rescaling 

in < n > Y + A 

where £ = hi(P/E) = ln(l/x) for a particle of energy E in a jet with primary parton mo- 
mentum P,Y = ln(PO/Qo) and A = ln(Qo/A). The function f(() has an approximately 
Gaussian shape, the so-called "hump-backed plateau" PJl9|. Again, the approach to this 
limit is rather slow, for example, the maximum of the spectrum occurs at (for Qq « A) 

r = n^ + ^-^) + o.i (io) 

where c = 0.2915 for rif = 3. The leading DLA term gives the asymptotic limit (* — > |, 
the next two terms the high energy corrections | 20| , the last term a numerical estimate 



of the remaining contributions applicable in the present energy range jnj]. According 

to the LPHD hypothesis this prediction at the parton level can be compared directly to 

the hadronic observable. The rescaled quantity (* = ^* /Y for charged hadrons is shown 

in Fig. 1 up to LEP-1.5 energies in comparison with the prediction ([H]). The cut-off 

parameter Qq = 0.270 GeV is taken from a global fit to the moments of the distribution 

from cms energies 3 to 91 GeV in e + e~ annihilation [p2| . The data in Fig. 1 closely 

follow the MLLA prediction ( |10"D which very slowly approaches the asymptotic DLA limit 
/•* i 

^ ~ 2' 



3.3 Angular correlations 



There was a lot of interest in the last years in the study and interpretation of angular 
correlations | 24| , which was triggered by the suggestion p5[ , such correlations could be 
power behaved at high resolution ( "intermittency" ) . Such power behaviour is expected, 
for example, for selfsimilar branching processes, so it applies to QCD to the extent that 



arc 



the running of the coupling is neglected f2"8| . 

The observables which are considered in the analytical QCD calculations P2] , pB|P? 
the two particle correlation density between two particles p ( - 2 - ) ($i2, O) in the forward cone 
of half angle 0, the factorial or cumulant multiplicity moments and for particles 
between two cones at angles — 5 and Q + 5 around the jet axis or in a cone of angular size 
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6 in direction with respect to the jet axis. According to the volume of phase space 5 D 
these two configurations are refered to by their dimensions D = 1 und D = 2. One finds 
the results on these correlations by first deriving the integral equation of the n-particle 
correlation function from @ and @, and then integrates over the remaining variables; 
the resulting evolution equation is of the type ([]) which can be solved approximately for 
running a s (for fixed a s one can get often exact results). 

The correlation functions of order n are conventionally normalized by a power of the 



multiplicity < n > ~ exp(2/3y ln(P0/A)). These normalized correlations, after removal 
of certain known kinematic phase space factors follow, after rescaling, the asymptotic 



angular scaling law 12 



ln#( n )(5,9,P) In 



2n(3J\nPe/A 



(l-u(e,n)/n), e = (11) 

111 A 



in the rescaled logarithmic variable e with < £ < 1. So the rescaled observables do 
not depend on the variables 5, G and P separately but only through the variable e. The 
normalization of the l.h.s of ( p"T|) corresponds to (In < n >) n . The function u; is known 
analytically, for large n one finds 

u(e,n) = n v / T^I(l-^ ln(l - e) + . . .) (12) 

which turns out to be a good approximation already for n — 2. 

An interesting feature of these results is their universality, i.e. the same limit is ob- 
tained for quite different observables: the correlations f($i 2 ) = p^ 2 -* (^12, @) / < n(O) > 2 
and the normalized moments of any order in one or two dimensions. These correlation 
functions refer actually to particles in quite different regions of phase space. 

As an example we show in Fig. 2 the rescaled normalized two particle density f ($12, 6) 



as obtained from the DELPHI collaboration [29| which is rather well approximated by the 



asymptotic prediction from the DLA. The data are in good agreement with the Monte 
Carlo calculations at the same energy at either parton or hadron level. Monte Carlo 
results at a much higher energy show a similar behaviour in agreement with the scaling 
prediction (|TTp . An observable which projects out the genuine 2-particle correlations more 
effectively from the uncorrelated background is the "correlation integral" [30,12] r($i2) = 
p^ 2 \'d 12) / Pnorm($ '12) where the normalization corresponds to the density of relative angles 
$12 of particles from different jets. This quantity has been measured as well and the 



predicted angular scaling law (11]) for r^^) has been verified for different jet opening 
angles 0. 

An uncertainty in these comparisons, which is hard to quantify, comes from the choice 
of the jet axis (taken usually as the sphericity axis). An improvement of both theoretical 
and experimental results could be obtained by using the Energy-Multiplicity-Multiplicity 
(EMM) definition as applied already to the 2-particle azimuthal angle correlations pi| . 



The predictions for moments have been derived for cumulants |L2| or for factorial 
moments pB|P?| which approach the same limit asymptotically. It turns out that the 
factorial moments at present energies are much closer to the asymptotic predictions; 
such moments are shown in Fig. 3 after appropriate rescaling by kinematic factors again 
as function of the e-variable. The data show the same trend as the asymptotic DLA 
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predictions of fllTI). The approach of the cumulant moments to the asymptotic results is 
much slower. 

Note that the curve in Fig. 3 for n = 2 is the same as the one in Fig. 2 according to the 
universality property of (0). Also the differences between the moments of orders n = 2 
and n = 3 are largely removed after rescaling; they are of relative order 1/n 2 according 
to 



11) and 01 



The various angular regions have quite different characteristics. For small e (large 
relative angles) the function u(e, n) « n — \ IL -^-£, which yields a power behaviour of the 
moments (either C {n) or F^) 



AfW(^5)~f|j , n = D(n-l)-(n-i) T o(P0). (13) 

This is the asymptotic power law ( "intermittency" ) for the QCD cascade. It applies for 
large relative angles where the cascade is fully developed and reflects the selfsimilarity of 
the branching process; the intermittency exponent <fr n depends on the scale through the 
running a s . In case of fixed a s the same result is obtained with universal 70 parameter. 
Moving to larger e the observables show the angular scaling law (|TT|) with non-linear 
function u(e,n). In this region the results are infrared safe (i.e. do not depend on Qo), 
they apparently are also not much dependent on hadronization effects (see Fig. 2). 

Moving to yet larger e one comes to a critical angle e cr u which separates two kinematic 



regimes of quite different characteristics [28[31[32| . The correlation functions have a 



discontinuous second derivative at this angle. In the new region e > e cr i t (small relative 
angles) the correlation functions do depend on the cut-off Q and they become independent 
of the order n, contrary to the behaviour for e < e cr u] they are given in terms of the one- 
particle inclusive spectrum. As a consequence, one expects in this region a dependence 
of the particle type (titi, KK, pp correlations) if the particle mass is related to the cut-off 

Qo- 



For fixed a s this angle is given by e cr i t = n /(n + 1) at order n |pl| , p2| . For running 



a s a new scale A appears and one considers p2[ the double scaling limit of the function 



oo n (e } p) with p = JX/(Y + A) for asymptotic Y at fixed e and p. Then again a critical 
behaviour at a certain angle e cr n is found. This limit requires also an increasing Q 
and therefore does not correspond to the usual fixed cut-off. Therefore, for finite, 
physical Q the separation of the two regions is not expected to be complete. It would be 
interesting to verify the characteristics of these two regimes experimentally. 



4 Scaling law for soft particles 

The prediction of the hump-backed shape of the inclusive energy spectrum in the £ = 
lnl/x variable and its subsequent observation was an important success of QCD in its 
application to multiparticle physics. The coherence of the soft gluon emission from all 
harder partons forbids the multiplication of the soft particles and one expects nearly an 
energy independence of the soft particle rate [[| . Such a property has been pointed out 
to be present indeed in the data [^ , p2]| . 

This problem has been studied recently in more detail [3TJ. The analytical calcula- 
tions both in DLA and MLLA converge towards the same limits independent of the cms 
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energy for small particle energies. In this limit the energy conservation effects and large 
z corrections from the splitting functions which make up the differences between the ap- 
proximations (|3|) and (H) can be neglected. If LPHD is valid towards these low energies 
one expects also a scaling behaviour for the invariant density Iq of hadrons in the soft 
limit where the particle momentum p or rapidity y and transverse momentum fc^ become 
small: 

I = hm E— = -limE— . (14) 
y^o,p T ^o d A p 2 p^o d s p 

The factor | in this definition takes into account that both hemispheres are included in the 
limit p — > 0. This scaling behaviour is a direct consequence of the coherence of the gluon 
emission: The emission rate for the gluon of large wavelength does not depend on the 
details of the jet evolution at smaller distances; it is essentially determined by the colour 
charge of the hard initial partons and is energy independent. The energy independent 
contribution comes from the single gluon bremsstrahlung of order a s , the higher order 
contributions generate the energy dependence but do not contribute in the soft limit. 

In Fig. 4 we show the experimental results on the invariant density of charged parti- 
cles for cms energies from 3 to 130 GeV in e + e~ annihilation. An approximate energy 
independence of the soft limit (within about 20%) is indeed observed; the same is true 
for identified particles it, K and p ||34j| . The curves in Fig. 4 represent the MLLA results, 
where also a particular prescription is employed to relate the different parton and hadron 
kinematics near the boundary E ~ Qo (for more details, see |34| , |3"5|| ). The theoretical 
curves show the approach to the scaling limit and describe well the different slopes at 
larger particle energies. An important role here is played by the running a s which pro- 
vides the strong rise towards small energies for E < 1 GeV, for fixed a s this rise would 
be much weaker p2| , |34| . 

A crucial test of the QCD-LPHD interpretation of this scaling result is the verification 
of the dependence of the limiting densitiy I on the primary colour charge. This can 
be obtained from e + e~ — >3 jets, deep inelastic scattering or semihard hadronic processes 



with gluon exchange [34 



5 Summary 

The perturbative approach to multiparticle production in connection with the LPHD 
assumption represents a very economic description of the phenomena which involves only 
the parameters Qo and A apart from the normalization. 

The analytical treatment singles out the logarithmic momentum and angular variables 
which are appropriate to the description of bremsstrahlung processes and absorb the 
collinear and soft divergent behaviour. Therefore the finite asymptotic limits of various 
observables in the rescaled logarithmic variables are a direct consequence of the parton 
branching process generated by bremsstrahlung type emissions. These scaling laws are 
then more specific to QCD than the KNO multiplicity scaling which holds for a wide 
class of branching processes, not necessarily of bremsstrahlung type. These results are 
obtained in the DLA where energy conservation is neglected. A noteworthy feature of 
angular correlations not met in energy spectra is the occurence of a critical angle which 
separates two scaling regimes with quite different characteristics. 
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Another scaling prediction from DLA is obtained in the soft particle limit at finite 
energies where energy conservation effects can be neglected as well. It is remarkable that 
perturbative QCD predictions work even in such an extreme limit and this requires further 
investigations with different partonic antenna patterns for confirmation. 



Acknowledgements 

I would like to thank A. Bialas for his inspiration of the studies of scaling laws in mul- 
tiparticle physics and V. A. Khoze, S. Lupia and J. Wosiek for the collaboration on the 
subjects of this lecture. 



8 



References 

[I] T. Sjostrand, Comp. Phys. Commun., 39 (1986) 347; 

T. Sjostrand and M. Bengtsson, Comp. Phys. Commun., 43 (1987) 367. 

[2] G. Marchesini and B. R. Webber, Nucl. Phys., B238 (1984) 1; 310 (1988) 461; 
G. Marchesini, B. R. Webber, G. Abbiendi, I. G. Knowles, M. H. Seymour, and L. 
Stanco, Comp. Phys. Comm. 67 (1992) 465. 

[3] Ya. I. Azimov, Yu. L. Dokshitzer, V. A. Khoze and S. I. Troyan, Z. Phys., C27 (1985) 
65 and C31 (1986) 213. 

[4] Yu. L. Dokshitzer, V. S. Fadin and V. A. Khoze, Phys. Lett., 115B (1982) 242; 
Z. Phys., C15 (1982) 325. 

[5] A. Bassetto, M. Ciafaloni and G. Marchesini, Phys. Rep., C100 (1983) 201. 

[6] A. H. Mueller, Nucl. Phys., B213 (1983) 85; Erratum quoted ibid., B241 (1984) 141; 
Yu. L. Dokshitzer and S. I. Troyan, Proc. 19th Winter School of the LNPI, Vol. 1, 
p. 144; Leningrad preprint LNPI-922 (1984). 

[7] Yu. L. Dokshitzer, V. A. Khoze, A. H. Mueller and S. I. Troyan, "Basics of Perturbative 
QCD", ed. J. Tran Thanh Van, Editions Frontieres, Gif-sur-Yvette, 1991. 

[8] for a recent discussion of the phenomenological status of the analytical perturbative 
approach, see V. A. Khoze and W. Ochs, 'Perturbative QCD approach to multiparticle 
production", preprint Durham DTP/96/36, MPI-PhT-96/29 (1996). 

[9] B. I. Ermolayev and V. S. Fadin, JETP. Lett., 33 (1981) 285; 
A. H. Mueller, Phys. Lett., 104B (1981) 161. 

[10] A. H. Mueller, Phys. Rev., D4 (1971) 150. 

[II] Yu. L. Dokshitzer, V. S. Fadin and V. A. Khoze, Z. Phys., C18 (1983) 37. 

[12] W. Ochs and J. Wosiek, Phys. Lett., B304 (1992) 144; Z. Phys., C68 (1995) 269. 

[13] Z. Koba, H.B. Nielsen and P. Olesen, Nucl. Phys., B40 (1972) 317. 

[14] A.M. Polyakov, Sov. Phys. JETP, 32 (1971) 296; 33 (1971) 850. 

[15] S. J. Orphanidis and V. Rittenberg, Phys. Rev., D10 (1974) 2892. 

[16] A. Bassetto, M. Ciafaloni, G. Marchesini, Nucl. Phys., B163 (1980) 477. 

[17] E. D. Malaza, B. R. Webber, Phys. Lett., 149B (1984) 501; Nucl. Phys., B267 
(1986) 70. 

[18] W. Ochs, Invited talk at the XXIV Int. Symp. on Multiparticle Dynamics, Vietri 
sul Mare, Italy, Sept. 1994, eds. A. Giovannini, S. Lupia, and R. Ugoccioni, (World 
Scientific, Singapore), p. 243. 



9 



[19] A. Bassetto, M. Ciafaloni, G. Marchesini and A. H. Mueller, Nucl. Phys., B207 
(1982) 189. 



Yu. L. Dokshitzer, V. A. Khoze and S. I. Troyan, Int. J. Mod. Phys., A7 (1992) 1875. 
[21] V. A. Khoze, S. Lupia and W. Ochs, Phys. Lett., B386 (1996) 451. 
[22] S. Lupia and W. Ochs, Phys. Lett., B365 (1996) 339. 
[23] S. Lupia and W. Ochs, in preparation. 

[24] E. A. De Wolf, I. M. Dremin and W. K. Kittel, Phys. Rep. 270 (1996) 1. 

[25] A. Bialas and R. Peschanski, Nucl. Phys., B273 (1986) 703; B308 (1988) 857. 

[26] Yu. L. Dokshitzer and I. Dremin, Nucl. Phys., B402 (1993) 139. 

[27] Ph. Brax, J. L. Meunier and R. Peschanski, Z. Phys., C62 (1994) 649. 

[28] W. Ochs and J. Wosiek, Phys. Lett., B289 (1992) 159. 

[29] DELPHI Collaboration (preliminary results): B. Buschbeck, F. Mandl et al., Con- 
tributed paper EPS 0553 to 1995 Int. Europhysics Conference on High Energy Physics, 
27 July - 2 August 1995, Brussels, Belgium, DELPHI 95-97, PHYS 532 (1995). 

[30] P. Lipa, P. Carruthers, H.C. Eggers and B. Buschbeck, Phys. Lett., B285 (1992) 
300. 

[31] Yu. L. Dokshitzer, G. Marchesini and G. Oriani, Nucl. Phys., B387 (1992) 675. 
[32] W. Ochs and J. Wosiek, Z. Phys., C72 (1996) 263. 

[33] V. A. Khoze, in in "Z° Physics", Cargese 1990, Eds. M. Levy, J-L. Basdevant, M. 
Jacob, D. Speiser, J. Weyers and R. Gastmans, NATO ASI, Series B: Physics Vol. 261, 
p. 419. 

[34] V. A. Khoze, S. Lupia and W. Ochs, "QCD coherence and the soft limit of the energy 
spectrum", preprint DTP/96/90, MPI-PhT/96-92, |iep-ph/9610"20l October 1996. 

[35] V. A. Khoze, S. Lupia and W. Ochs, "Soft particle production and QCD coherence", 
preprint DTP/96/94, MPI-PhT/96-94, |hep-ph/961034^ , October 1996. 



10 



1.0 



0.8 - 



C 



* 0.6 




0.4 



0.2 



0.0 



♦ TASSO 
ill MARK II 
+ TPC 

* CELLO 
X AMy 



-M TOPAZ 
® ALEPH 
□D DELPHI 
X L3 
CD OPAL 



5 

Y 



6 



Figure 1: Maximum of the rescaled inclusive momentum distribution (* = (*/^ as a 
function of Y — hi 7^; comparison between experimental data from e + e~ annihilation 
and theoretical prediction in MLLA numerically extracted from the shape of the Limiting 
Spectrum (solid line) for the cut-off parameter Q = A = 270 MeV. Crosses mark the 
predictions at the cms energies 200 GeV and 500 GeV. Asymptotically, the leading DLA 
result (* = I is approached (see plU2~3f ). 
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Figure 2: The rescaled 2-particle angular correlation f = p^i'&n)/ < n > in the 
forward cone with half-opening O as function of the scaling variable e as measured by 
DELPHI [g9|. Also shown are the results from the JETSET and HERWIG Monte Carlo's 
at the parton and hadron levels at different energies. The data show the predicted scaling 
behaviour and the approach to the asymptotic DLA prediction (for A = 0.15 GeV, with 

b = 2/3 v /ln(P6/A)). 
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Figure 3: Rescaled factorial multiplicity moments for particles in the ring around the jet 
axis with polar angles between — $ and + i9 as measured by DELPHI |29| (momentum 
P = 45 GeV, A = 0.15 GeV) in comparison with the asymptotic DLA prediction. Note 
that the curve for n = 2 is the same as in Fig. 2 for the rescaled correlation f. 
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Figure 4: Invariant density Edn/d 3 p of charged particles in e + e annihilation as a function 

of the particle energy E = \Jp 2 + Qq at Q = 270 MeV. Data at various cms energies are 
compared to MLLA predictions with the overall normalization adjusted (from [j34T]). 
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